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ABSTRACT 



Aims. To reexamine the implications of the recent HESS observations of the blazar 1ES0229+200 for constraining the 
extragalactic mid-infrared background radiation. 

Methods. We examine the effect of 7-ray absorption by the extragalactic infrared radiation on predicted intrinsic spectra 
for this blazar and compare our results with the observational data. 

Results. We find agreement with our previous results on the shape of the IR spectral energy distribution (SED), contrary 
to the recent assertion of the HESS group. Our analysis indicates that 1ES0229+200 has a very hard intrinsic spectrum 
with a spectral index between 1.1 ± 0.3 and 1.5 ± 0.3 in the energy range between ~0.5 TeV and ~15 TeV. 
Conclusions. Under the assumptions that (1) the SED models of Stecker, Malkan & Scully (2006) are reasonable as 
derived from numerous detailed IR observations, and (2) spectral indexes in the range 1 < V < 1.5 have been shown 
to be obtainable from relativistic shock acceleration under the astrophysical conditions extant in blazar flares (Stecker, 
Baring & Summerlin 2007), the fits to the observations of 1ES0229+200 using our previous IR SEDs are consistent with 
both the IR and 7-ray observations. Our analysis presents evidence indicating that the energy spectrum of relativistic 
particles in 1ES0229+200 is produced by relativistic shock acceleration, producing an intrinsic 7-ray spectrum with 
index 1 < T < 1.5 and with no evidence of a peak in the SED up to energies ~ 15 TeV. 
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1. Introduction 

Shortly after the first strong 7-ray blazar 3C279 was 
discovered by the EGRET detector aboard the Compton 
Gamma Ray Observatory (CGRO), Stecker, De Jager & 
Salamon (1992) proposed that the study of the spectra of 
such sources could be used to probe the intergalactic in- 
frared radiation. At that time, there were no actual direct 
observations of the diffuse extragalactic IR background or 
extensive observations of the sources of such radiation. The 
idea was to look for the effects of photon-photon annihi- 
lation interactions into electron-positron pairs. The cross 
section for this process is exactly determined; it can be cal- 
culated using quantum electrodynamics (Breit & Wheeler 
1934). Thus, in principle, if one knows the emission spec- 
trum of an extragalactic source at a given redshift, one 
can determine the column density of photons between the 
source and the Earth. 

In the last 15 years, great advances have been made in 
extragalactic infrared astronomy. The diffuse background 
at wavelengths not totally dominated by galactic or zodia- 
cal emission has been measured by the Cosmic Background 
Explorer (COBE). In addition, there have been extensive 
observations of infrared emission from galaxies themselves, 
whose total emission is thought to make up the cosmic in- 
frared background (see review by Hauser & Dwek 2001). 
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The latest extensive observations have been made by the 
Spitzer satellite. It is thus appropriate to use a synoptic 
approach combining the TeV 7-ray observations with the 
extragalactic infrared observations in order to best explore 
both the TeV emission from blazars and the diffuse extra- 
galactic infrared radiation. 

Aharonian et al. (2007) have recently observed the spec- 
trum of the BLLac object 1ES0229+200 up to an energy 
~ 15 TeV with the High Energy Spectroscopic System 
(HESS). Then, by assuming that the intrinsic spectral in- 
dex of this source is greater than 1.5, they drew conclusions 
regarding wavelength dependence and flux of the mid-IR 
extragalactic background radiation. Their conclusions re- 
garding the mid-IR extragalactic background radiation ap- 
pear to disfavor the results of the extensive semi-empirical 
calculations of the extragalactic IR background spectrum 
given by Stecker, Malkan & Scully (2007) (SMS). 

In this paper, we will reexamine the assumptions and 
conclusions of Aharonian et al. (2007) and show that the 
observations of 1ES0229 are fully consistent with the diffuse 
IR background spectrum obtained by SMS. Furthermore, 
we will show that 1ES0229+200 is an example of a set of 
blazars which exhibit very hard spectra that are indicative 
of relativistic shock acceleration. 
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2. The Diffuse Extragalactic Infrared Background 

Various calculations of the extragalactic infrared back- 
ground have been made (Stecker, Puget & Fazio 1977; 
Malkan & Stecker 1998, 2001; Totani & Takcuchi 2002; 
Kneiske, et al. 2004; Primack, Bullock & Somervillc 2005; 
Stecker, Malkan & Scully 2006, hereafter designated SMS). 
Of these models, the most empirically based are those of 
Malkan & Stecker (1998,2001), Totanti & Takeuchi (2002) 
and SMS. The SMS calculation includes input from the lat- 
est Spitzer observations. Because the largest uncertainty in 
these calculations arises from the the uncertainty in the 
temporal evolution of the star formation rate in galax- 
ies, SMS assumed two different evolution models, viz., a 
"baseline" model and a "fast evolution" model. These mod- 
els produced similar wavelength dependences for the spec- 
tral energy distribution (SED) of the extragalactic infrared 
background, but give a difference of roughly 30-40% in over- 
all intensity. 

The empirically based calculations mentioned above in- 
clude the observationally based contributions of warm dust 
and emission bands from polycyclic aromatic hydrocarbon 
(PAH) molecules and silicates that have been observed to 
contribute significantly to galaxy emission in the mid-IR 
(e.g., Lagache et al. 2004). These components of galac- 
tic infrared emission have the effect of partially filling in 
the "valley" in the mid-IR SED between the peak from 
starlight emission and that from dust emission. In contrast, 
the model of Primack et al. (2005) exhibits a steep mid-IR 
valley, which is in direct conflict with lower limits obtained 
from galaxy counts at 15 /xm (Altieri et al. 1999; Elbaz et 
al. 2002). 



3. Intrinsic Spectra of Blazars from Theoretical 
Considerations 

The key difference between our analysis and the analysis of 
Aharonian et al. (2007) is their assumption that the intrin- 
sic spectral indexes of blazars cannot be greater than 1.5. 
Stecker, Baring and Summcrlin (2007) (SBS) have shown 
that spectral indexes between 1 and 1.5 are obtainable from 
relativistic shock acceleration. They list three blazars with 
intrinsic spectral indexes between ~ 1 and ~ 1.5 in the en- 
ergy range between 0.2 TeV and 2 TeV. Their simulations 
indicate that a range of spectral indexes can be produced 
by relatvistic shocks, depending on the conditions in the 
individual shocks. 

In fact, recent observations of an extreme MeV 7-ray 
blazar at a redshift of ^3 by Swift (Sambruna et al. 2006) 
and the powerful 7-ray quasar PKS 1510-089 at a redshift of 
0.361 by Swift and Suzaku (Kataoka et al. 2007) both exhib- 
ited power-law spectra in the hard X-ray range that had in- 
dexes less than 1.5, implying electron spectra with indexes 
less than the value of 2 usually considered for shock acceler- 
ation. In addition, the quasar IGR J22517+2218, at a red- 
shift of 3.668, has been observed with IBIS/INTEGRAL to 
have a spectral index of 1.4 ± 0.6 in the 20-100 keV hard X- 
ray energy range (Bassani et al. 2007). Spectra in the hard 
X-ray range do not suffer intergalactic absorption so that 
there is no ambiguity concerning their spectral indexes. 



4. The Observed Spectrum of 1ES0229+200 and 
its Derived Intrinsic Spectrum 

The BL Lac object 1ES0229+200 was predicted to be an 
observable TeV source with a TeV flux of ~ 10~ 12 cm _2 s _1 
by Stecker, de Jager & Salamon (1996). It was detected by 
Aharonian et al. (2007) at a flux level of - 0.94 x 10" 12 
above 0.58 TeV. Aharonian et al. (2007) have measured 
the spectrum of the blazar 1ES0229+200 at energies be- 
tween 0.5 TeV and ~15 TeV. This blazar lies at a redshift 
of 0.14 (Schachter et al. 1993). This redshift, together with 
the spectral observations, implies a very hard intrinsic spec- 
trum for the source. 

According to the observations of Aharonian et al. 
(2007), the blazar 1ES0229+200 has an observed spectral 
index of 2.50. In their analysis, Aharonian et al. chose to 
"deabsorb" their data points using the results of various op- 
tical depth calculations and to fit a simple power-law spec- 
trum to the results. Using estimated optical depths from 
SMS, they gave best-fit power-law spectral indexes of 0.1 
± 0.3 for the fast evolution case and 0.6 ± 0.3 for the base- 
line case. We have done a reanalysis using their method and 
find spectral indexes of 0.2 ± 0.3 and 0.8 ± 0.3 for the SMS 
fast evolution and baseline models respectively. The small 
difference in the results can be attributed to Aharonian et 
al. using interpolated optical depths, as SMS did not pro- 
vide the results at the specific redshift of 1ES0229+200. 
We here calculated the optical depths for z = 0.14 exactly, 
for both the fast evolution and baseline evolution models 
of SMS. 

In this paper, we will adopt a different method for 
analysing the intrinsic spectrum of 1ES0229+200, one 
which we feel is superior to the above approach. Aharonian 
et al.(2007) chose to force a power-law fit to the deabsorbed 
data points. We choose here to assume an intrinsic power- 
law spectrum emitted by the source and then apply a cor- 
rection for optical depth. We then compare the resulting 
spectrum with the actual observed spectrum. Because of 
the nonlinear nature of the energy dependence of the op- 
tical depth, we do not expect that the observed spectrum 
will have a power-law form. 

Using the formulas given by Stecker & Scully (2006), we 
predict an intrinsic spectral index for 1ES0229+200 of 1.11 
for the fast evolution model SED and 1.45 for the baseline 
model SED in the 0.2 to 2 TeV energy range. This intrin- 
sic power-law spectrum, derived for the 0.2 TeV to 2 TeV 
energy range, is then assumed to be extended to ~15 TeV. 
Such a spectrum with an index between 1 and 1.5 can be 
obtained from relativistic shock acceleration, as discussed 
in the previous section. We have then calculated the form 
of the resulting spectrum that would be observed at Earth 
by taking account of the absorption predicted by both the 
baseline and fast evolution models of SMS. 

Thus, the observed spectrum is fixed to be of the form 
K FE E- 1A1 e- T ^ E ' z = ^ &ndK B E- 1A5 e- T ^ E ' z = QA ^ for 
the fast evolution and baseline models respectively, where r 
is the optical depth of the universe to 7-rays originating at 
a redshift of 0.14 calculated according to SMS. We employ 
a Levenberg-Marquardt nonlinear least squares method to 
fit this form of the spectrum to the observed data. In 
each case, the single parameter of the fit is the normaliza- 
tion coefficient, K. The Levenberg-Marquardt method min- 
imizes x 2 to best fit of a nonlinear function to the observa- 
tional data. For the fast evolution model, Kfe = 10~ 10 - 80 
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Fig. 1. The observed and calculated spectrum of 
1ES0029+200. Data are from Aharonian et al. (2007). The 
vertical bars are statistical errors; the horizontal bars are 
the energy bins. The calculations of the spectra are as de- 
scribed in the text. The solid line is the best fit with r given 
for the baseline model of SMS; the dashed line is the best 
fit with t given for the fast evolution model of SMS. 

cm _2 s —1 TeV _1 with E in TeV; for the baseline model, 
Kb = 10~ 1113 in the same units. The resulting spectra 
are plotted along with the data in Figure 1. The solid line 
is the best fit with r given for the baseline model of SMS; 
the steeper dashed line is the best fit with r given for the 
fast evolution model of SMS which gives larger values of r. 
As can be seen in Figure 1, these curves show excellent con- 
sistancy with the observational data obtained by the HESS 
group for 1ES0229 + 200. Indeed, the x 2 vlaues obtained 
are 3.3 for the baseline fit and 5. 7 for the fast evolution fit 
with 7 degrees of freedom. 

5. Conclusions and Discussion 

Aharonian et al. (2007), by assuming that blazar spectra 
have spectral indexes T > 1.5, concluded that the mid- 
IR spectral energy distribution (SED) must have a wave- 
length dependence steeper than A -1 , in their terminol- 
ogy A~ a ,a > 1.1 ± 0.25 in the wavelength range between 
2/im and 10/im. On the contrary, galaxy emission mod- 
els which take into account emission by warm dust and 
PAH and silicate emission in the mid-IR, as well as direct 
mid-IR observations of galaxy spectra give values for a in 
the range between ~ 0.7 and ~ 0.8 (Spinoglio et al 1995, 
Xu et al. 2001). These values are consistent with the 2- 
10/im diffuse background SED given in the semi-empirical 
SMS models. Such values lead to an energy dependence for 
the 7-ray optical depth which is consistent with that ob- 
tained by both Totani & Takeuchi (2002) and SMS. This 
result is also supported by the lower limit on the 15/im dif- 
fuse background flux obtained by galaxy counts of 3.3 ±1.3 
nW m _2 sr _1 (Altieri et al. 1999) which, under the conser- 
vative assumption that 80% of the mid-IR flux is resolved 
out (SMS), yields a value for the total diffuse background 
flux at 15/J,m of 4.1 ± 1.6 nW m _2 sr , higher than the 
upper limit of 3.1 nW m sr derived by Aharonian et al. 
(2007) under the assumption that T > 1.5. 

Under the assumptions that (1) the SED models of 
Stecker, Malkan & Scully (2006) (SMS) are reasonable as 
derived from numerous detailed IR observations, and (2) 
spectral indexes in the range 1 < T < 1.5 have been shown 



by Stecker, Baring & Scully (2007) (SBS) to be obtainable 
from relativistic shock acceleration under the astrophysical 
conditions extant in blazar flares, the fits to the HESS TeV 
observations of 1ES0229+200 using the SMS infrared SEDs 
are consistent with both the IR and 7-ray observations. 

The SBS simulations indicate with specific test runs 
that electron spectra with asymptotic spectral indexes be- 
tween 1.26 and 1.62 can be obtained from acceleration by 
relativistic shocks with bulk Lorentz factors between 10 and 
30 and these electrons can then Compton scatter to produce 
7-ray spectra with indexes ~1.1 and ^1.3. The simulations 
show that larger bulk Lorentz factors lead to flatter spectra. 
Such Lorentz factors of 50 or more have been implied from 
studies of specific flares in BL Lac objects (Konopelko et al. 
2003; Bagelman, Fabian & Rees 2007) so that the existence 
of highly relativistic shocks in such sources with indexes as 
flat as ~1 is not unreasonable. 

For the SMS fast evolution SED, the 7-ray index ob- 
tained is 1.11 and for the baseline SED, the index ob- 
tained is 1.45. Our analysis thus presents evidence for rel- 
ativistic shock acceleration in 1ES0229±200 that results 
in a very hard intrinsic 7-ray spectrum with no evidence 
of a peak in the 7-ray SED up to energies ~ 10 TeV. 
Unfortunately, there are no simultaneous observations at 
other wavelengths that can be used to model the flare that 
occurred in 1ES0229+200 and, in any case, TeV orphan 
flares have been observed in other BL Lac objects. However, 
we note that hard X-ray spectra have been seen in other 
sources as discussed in Section 3. We also note that SBS 
have derived spectra for three other BL Lac objects at red- 
shifts between 0.18 and 0.19 with indexes between 1 and 
1.5, viz., 1ES1218±30, 1ES1101-232, and 1ES0347-121. It 
therefore appears that a whole class of blazars or blazar 
flares may exhibit the hard-spectrum characteristics of rel- 
ativistic shock acceleration. 

Acknowledgements. We wish to thank Wystan Benbow for sending 
us a list of data on the spectrum of 1ES0229+200 observed by HESS. 
STS gratefully acknowledges partial support from the Thomas F. & 
Kate Miller Jeffress Memorial Trust grant no. J-805. We thank an 
anonymous referee for constructive comments. 



References 

Aharonian, F. et al. 2007, A&A, in press, larXiv:0709.4584l 
Altieri, B. et al. 1999, A&A 343, L65 

Bassani, L. et al. 2007, ApJ Letters, in press. larXiv70 709.3023 

Begelman, M. C, Fabian, A. C. & Rees, M. J. 2006, arXiv:0709.0540 

Breit, G. & Wheeler, J. A. 1934, Phys. Rev. 46, 1087 

Elbaz, D. et al. 2002, A&A 384, 848 

Hauser, M.G. & Dwek, E. 2001 ARA &A 39, 249 

Kataoka, J. et al. 2007, ApJ in press. larXTv:0709. 15281 

Kneiske, T. M. et al. 2004, A&A 413, 807 

Konopelko, A. et al. 2003, ApJ 597, 851 

Lagache,G. et al. 2004, ApJ Suppl. 154, 112 

Malkan, M. A. & Stecker, F. W. 1998, ApJ 496, 13 

Malkan, M. A. & Stecker, F. W. 2001, ApJ 555, 641 

Primack, J. R., Bullock, J. S. & Somerville, R. S. 2005, AIP Conf. 

Proc. 745, 23 
Sambruna, R. M. et al. 2006, ApJ 646, 23 
Schachter, J. F. et al. 1993, ApJ 412, 541 
Spinoglio, L. et al. 1995, ApJ 453, 616 

Stecker, F. W., Baring, M.G., & Summerlin, E.J. 2007, ApJ Letters 
667, L29 

Stecker, F. W., de Jager, O. C. & Salamon, M. H. 1992, ApJ Letters 
390, L49 

Stecker, F. W., de Jager, O. C. & Salamon, M. H. 1996, ApJ Letters 
473, L75 

Stecker, F. W., Malkan, M.A., & Scully S.T. 2006, ApJ, 648, 774 



4 



F.W. Stecker and S.T. Scully: 1ES0229 Spectrum and IR Background 



Stecker, F. W., Puget, J.-L. & Fazio, G. G. 1977, ApJ Letters 214, 
L51 

Stecker, F W. & Scully, S. T. 2006, ApJ Letters 652, L9 
Totani, T. & Takeuchi, T. T. 2002, ApJ 570, 470 
Xu, C. ct al. 2001, ApJ 562, 179 



